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Planetary and stellar dynamos likely result from turbulent motions in magnetofluids with kine¬ 
matic viscosities that are small compared to their magnetic diffusivities. Laboratory experiments are 
in progress to produce similar dynamos in liquid metals. This work reviews recent computations of 
thresholds in critical magnetic Reynolds number above which dynamo amplification can be expected 
for mechanically-forced turbulence (helical and non-helical, short wavelength and long wavelength) 
as a function of the magnetic Prandtl number Pm- New results for helical forcing are discussed, 
for which a dynamo is obtained at Pm = 5 x 1CU 3 . The fact that the kinetic turbulent spectrum 
is much broader in wavenumber space than the magnetic spectrum leads to numerical difficulties 
which are bridged by a combination of overlapping direct numerical simulations and subgrid models 
of magnetohydrodynamic turbulence. Typically, the critical magnetic Reynolds number increases 
steeply as the magnetic Prandtl number decreases, and then reaches an asymptotic plateau at values 
of at most a few hundred. In the turbulent regime and for magnetic Reynolds numbers large enough, 
both small and large scale magnetic fields are excited. The interactions between different scales in 
the flow are also discussed. 

PACS numbers: 47.65.+a; 47.27.Gs; 95.30.Qd 


I. INTRODUCTION 

Plasmas in stellar interiors and conducting fluids in 
planetary cores are characterized by a magnetic Prandtl 
number Pm (the ratio of the kinematic viscosity v to 
the magnetic diffusivity rf) much smaller than one. As a 
few examples, the magnetic Prandtl number in the solar 
convective region is estimated to be Pm ~ 10 -5 — 10 -6 
Q, and in the Earth’s core Pm ~ 10 -5 . Liquid sodium 
experiments are also characterized by small values of Pm- 

While numerical simulations of dynamo action in these 
objects are available, the large values of the kinetic ( Ry ) 
and magnetic ( Rm ) Reynolds numbers forbid a study us¬ 
ing realistic values of Pm- Simulations of the geodynamo 
[3 or the solar convective region |3j are often done for 
Pm ~ 1. While the proper separation of the kinetic and 
magnetic dissipation scales cannot be achieved in these 
simulations, values of Pm much smaller than one can be 
reached under more idealized conditions. Pseudospec- 
tral methods in periodic boxes give an excellent tool to 
study the behavior of magnetohydrodynamic (MHD) tur¬ 
bulence in the regime Pm < 1- The further assumption 
of incompressibility allows for an extra gain in computer 
power. Paralellized pseudospectral codes have reached 
for hydrodynamic turbulence resolutions of 4096 3 grid 
points, and Taylor Reynolds numbers of R\ « 1200 Q. 
These methods are conservative and nondispersive, be¬ 
ing well suited for the exploration of turbulent flows in 
regimes hard to explore in the laboratory. And under 
some circumstances, the range of values of Pm can be 
extended using subgrid scale (SGS) models. 

In this work, we review recent results from simulations 
of helical and non-helical dynamos at Pm < 1 using pseu¬ 
dospectral codes in periodic boxes [3, j6j, 0|. To extend 
the range in Pm in the simulations, SGS models were 
used in these works. We discuss some of these mod¬ 


els with particular emphasis in the Lagrangian Average 
MHD (LAMHD) equations |x. !i. 101:11111211. In addition, 
and to validate results from SGS models, new results 
from direct numerical simulations (DNS) with resolutions 
of 1024 3 grid points are presented. 

We focus on the properties of the magnetic field in the 
kinematic dynamo regime, when the intensity of the mag¬ 
netic field is small and the effect of the Lorentz force can 
be neglected. In particular, we discuss the behavior of the 
threshold in Rm for dynamo action as Pm is decreased. 
That is to say, given a hydrodynamic state and an arbi¬ 
trary small magnetic perturbation, what is the minimum 
value of Rm (or maximum value of rj, in some convenient 
set of dimensionless units) to have a dynamo instability 
such that the system reaches after a finite time a magne¬ 
tohydrodynamic steady state. Below this threshold, the 
magnetic perturbation is damped and the final state of 
the system is hydrodynamic. 

Two flows have recently been studied in this context: 
the flow resulting from Taylor-Green forcing III , and 
the result of Roberts forcing Hill- The first case corre¬ 
sponds to a flow with no net helicity that gives large scale 
dynamo action, while the former studies are for a helical 
flow where only small scale dynamo action is permitted 
by introducing mechanical energy in the largest available 
scale. In a different context, for isotropic, homogeneous, 
and delta-correlated in time external forcing, the prob¬ 
lem has also been studied in Refs. EMI. In addition to 
reviewing this results, we compare them against new sim¬ 
ulations using Arn’old-Childress-Beltrami (ABC) forcing 
with energy injection at intermediate scales. This is a 
helical case where large scale magnetic amplification is 
allowed. For this forcing, values of Pm down to 5 x 10~ 3 
are reached. The results obtained for such a low value 
of Pm are expected to be of relevance for astrophysical 
and geophysical applications, as well as for laboratory 
dynamos. 
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For all cases where a large scale flow is present, dynamo 
action is observed to persist at the smallest values of P M 
that can be reached. Moreover, for values of Pm smaller 
than ~ 0.1 an independence of the threshold with Pm is 
observed. While for the Taylor-Green (non-helical) forc¬ 
ing and the Roberts forcing (helical, but with magnetic 
amplification only at small scales) a sharp increase in the 
critical parameter is observed before reaching the asymp¬ 
totic regime, in the ABC case almost no such increase is 
found. 

The structure of the paper is as follows. In Sec. QH we 
present the equations, and the several forcing functions 
used. We also describe the code and the SGS models. 
Section lnTl discusses the thresholds for dynamo action for 
the several flows, and Sec. E] discusses the role played 
by the different scales in the problem. Finally, Sec. El 
presents the conclusions. 

II. EQUATIONS 
A. The MHD equations 

In dimensionless Alfvenic units, the incompressible 
MHD equations are 

— + v-Vv = -VP+j x B + i/V 2 v + f, (1) 

HD 

—+v-VB = B Vv + r)V 2 B, (2) 

with V • v = V • B = 0. Here v is the velocity field, and 
B is the magnetic field, related to the electric current 
density j by j = V x B. V is the normalized pressure-to- 
density ratio, obtained by solving the Poisson equation 
that results from taking the divergence of Eq. m and 
using the condition V • v = 0. 

In Eq. 0 f is an external force. Since in the incom¬ 
pressible framework we lack the usual energy sources in 
astrophysics and geophysics (e.g. thermal convection), 
we will set f to generate a large scale flow, and let the 
instabilities of the flow generate turbulent fluctuations. 
The expressions for f considered are discussed in the fol¬ 
lowing subsection. 

Equations 0 and 0 are solved using a standard pseu- 
dospectral method IJ, 0 in a three dimensional peri¬ 
odic box. The code uses projection into the Fourier base 
to compute spatial derivatives, and Runge-Kutta of ad¬ 
justable order to evolve the equations in time. The 2/3- 
rule for dealiasing is used. As a result, if N grid points 
are used in each direction, the maximum wavenumber 
solved by the code is kmax = AT/3. To make use of par¬ 
allel computers, the three dimensional Fourier transform 
has to be parallelized efficiently, using a methodology as 
described e.g. in Refs. .18, 19]. 

All Reynolds numbers discussed are based on the flow 
integral scale 

L = 2tt [ E v {k)k~ l dk/E v (3) 


where Ey{k) is the kinetic energy spectrum, and Ey = 
J Ey{k)dk is the total kinetic energy. Given the r.m.s. 
velocity U, the kinetic Reynolds number is Ry = LU/v 
and the magnetic Reynolds number is Rm = LU/rj. The 
magnetic Prandtl number is Pm = Rm/Rv- We also 
define the large scale eddy turnover time as T = L/U. 
The kinetic and magnetic dissipation wavenumbers in 
the turbulent regime are given by k v = (e/i / 3 ) 1 / 4 and 
krj = (e/ 77 3 ) 1 / 4 respectively, where e is the energy injec¬ 
tion rate. In all simulations, these wavenumbers were 
smaller than the maximum resolved wavenumber k max . 

The strategy is to turn on the forcing at t = 0 and allow 
the code to run for a time as a hydrodynamic code, with 
the magnetic field set to zero. Once a stationary state 
is reached, the magnetic field is seeded with randomly- 
chosen Fourier coefficients and allowed to amplify or de¬ 
cay. 


B. External forcing 

Several expressions for the external forcing were stud¬ 
ied in the context of dynamo action. Refs. 0. 0 consid¬ 
ered Taylor-Green [20j forcing 

f T G = [sin(kTGx)cos(k TG y)cos(k TG z)x- 

- cos(k TG x) sin{k TG y) cos(k TG z)y} , (4) 

with wavenumber /cxg = 2 (this election gives a peak in 
the kinetic energy at k « 3). If Eq. 0 is used as an ini¬ 
tial condition for the velocity field v, it is found that the 
flow is not a solution of the Euler’s equation and though 
highly symmetric, it leads to the rapid development of 
small spatial scales. The forcing is also non-helical, in 
the sense that fxG • V x fxG is zero at every point in 
space. Helical fluctuations can however appear at scales 
smaller than the forcing scale. 

Ref. HO considered the Roberts flow mm as the 
expression for the external force f: 

fR = \g sin(fc R x) cos (k R y)x - g cos(fc R x) sin (k R y)y+ 

+2/ sin(fc R x) sin(fc R y)z], (5) 

with the choice /c R = 1 . Since in that case mechanical 
energy is injected in the largest available scale, magnetic 
excitations can only grow at scales smaller than the en¬ 
ergy injection scale. The coefficients / and g are arbitrary 
and their ratio determines the extent to which the flow 
excited will be helical. We will concentrate here upon 
the case / = g. In this case, the forcing injects helicity 
in the flow, although not maximally. If used as an initial 
condition, the velocity field is an exact solution of Euler’s 
equation. As a result, the development of small scale fluc¬ 
tuations only takes place above a certain threshold in the 
kinetic Reynolds number Ry. For values of Ry smaller 
than this threshold, there is an exact laminar solution of 
the hydrodynamic equations given by v = f/(i/fc^). 
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We compare the results from these two forcing func¬ 
tions with the ABC forcing 


fABC = {[-Bcos(fc A BC2/) + Csin(fc A BC-z)]£+ 

+ [A sm(k ABC x) + C cos(fc AB c-z)] V + 

+ [A cos(fc AB c*) + B sin(fc AB cy)] z} , (6) 

with k A bc = 3, and A = 0.9, B — 1, and C = 1.1. 
Previous results for Pm < 1 have been reported in Refs. 
0 0, although no systematic study of the effect of 
lowering Pm was done. 


incompressible case, the LAMHD equations are 


— + u s • Vv = -VjVttj, - VP + j x B s 

+^V 2 v + f, (7) 

+ u s ■ VB S = B s • Vu s + t?V 2 B.. (8) 


The pressure V is to be determined, as before, from the 
relevant Poisson equation. The energy in this system is 
given by E = /(v • u s + B ■ B s )/2 d 3 x. The subindex s 
denotes smoothed fields, related to the unsmoothed fields 
by 


The ABC flow is maximally helical, and the election of 
^abc allows for the growth of magnetic energy at scales 
larger than the energy injection scale. As for the Roberts 
flow, the ABC flow is also an exact solution of the Eu¬ 
ler’s equation, and turbulent fluctuations are generated 
as the result of a hydrodynamic instability at values of 
Rv larger than a certain threshold. The particular elec¬ 
tion of the parameters A, B, and C is done to break the 
symmetry of the flow and reduce the value of the hydro- 
dynamic threshold 00. 


C. Subgrid models of MHD turbulence 


v = (l - a 2 V 2 ) u s (9) 

B = (l — a 2 V 2 ) B s . (10) 

In the context of the dynamo at low Pm, some others 
SGS models have been used. Ref. [f| used a modified 
LES 0| where only the velocity field at scales smaller 
than the magnetic diffusion scale was modeled using a 
turbulent effective viscosity dependent on the wavenum¬ 
ber. A similar SGS model has been used in Ref. 0 to 
study dynamo action with delta-correlated in time ran¬ 
dom forcing, but with the expression of the effective vis¬ 
cosity based on the Smagorinsky-Lilly model. Ref. 0 
also used hyperviscosity, although it should be remarked 
that this method is known to give wrong growth rates for 
the magnetic energy in the anisotropic case 0 . 


As the value of Pm is lowered down, the separation be¬ 
tween the Ohmic and viscous dissipation scales increases. 
As a result, the velocity field has more small scale struc¬ 
ture than the magnetic field. This imposes a stringent 
limit on the smallest Pm that can be reached using DNS. 
At some point, the small scales in the flow cannot be 
solved anymore, and some form of SGS modeling is re¬ 
quired. 

While SGS models of hydrodynamic turbulence have 
a rich history, models for MHD flows are still in their 
infancy (see e.g. 0 )- One of the main difficulties in 
MHD is that hypotheses often made in hydrodynamics 
(e.g. locality of interactions in Fourier space) are not nec¬ 
essarily true for magnetofluids. The general expression 
of the MHD energy spectrum is not known. And several 
regimes can be expected according to whether the sys¬ 
tem is mechanically or magnetically forced, whether the 
fields are statistically aligned or not, etc. 

In this work we will focus on the LAMHD equations 
(or MHD a-model). The hydrodynamic a-model (see e.g 
@, 0 and references therein) was validated against simu¬ 
lations in Ref. ;2J] . It differs from large eddy simulations 
(LES) in that the filter is applied to the Lagrangian of 
the ideal fluid, and as a result Hamiltonian properties 
of the system are preserved. It was extended to ideal 
MHD [a 0L and tested against DNS of MHD turbu¬ 
lence pjilia. In the context of magnetoconvection, the 
LAMHD equations were also studied in Ref. 0. In the 


III. BEHAVIOR OF THE DYNAMO WITH P M 

Figure ^ (bottom panel) shows the critical magnetic 
Reynolds number R C M as a function of Ry for Taylor- 
Green forcing. Crosses connected with solid lines are 
obtained using DNS, while crosses connected with dotted 
lines are from LAMHD simulations (the largest resolution 
used by each method was of 512 3 grid points). An overlap 
of the two methods for three values of Ry was used to 
verify the SGS model. In Ref. @ higher values of Ry 
were reached for this flow using LES. 

Above the threshold, magnetic field excitations are am¬ 
plified exponentially until the system reaches a fully de¬ 
veloped MHD regime. Below the threshold, magnetic 
field perturbations are damped and the system reaches 
after long times a hydrodynamic regime. 

The threshold is obtained from simulations in the fol¬ 
lowing way: in the steady state of a hydrodynamic sim¬ 
ulation at a given value of Ry, a small magnetic field is 
introduced and several simulations varying Rm are car¬ 
ried. In each MHD simulation, the exponential growth 
or decay rate a = d log Em /dt of the magnetic energy is 
measured. The value of Rm for which a = 0 defines R C M - 
In practice, R C M is bounded between two values of Rm 
that give positive and negative a, and the actual value 
of Rm is obtained from a best fit to the a vs. Rm curve 
(see e.g. Refs. @, 11 ). This procedure gives the value of 
R c m with an error of the order of 10%. 
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FIG. 1: Integral lengthscale of the flow as a function of Rv 
(above), and critical magnetic Reynolds for dynamo action 
R c m as a function of Rv (below) for the Taylor-Green forcing. 
Crosses connected with solid lines are obtained from DNS, 
while crosses connected with dotted lines are obtained using 
the LAMHD equations. 


Figure n (top panel) also shows the flow integral scale 
as Ry is increased. This scale is measured in the steady 
state of the hydrodynamic run that preceeds the MHD 
simulations. At low Ry the flow is laminar and L is 
close to the forcing scale 27r/fcTG- However, as Ry is 
increased, more and more small scale excitations appear 
in the flow, and L decreases. For Ry « 1000, the system 
has reached a fully turbulent regime with a Kolmogorov’s 
energy spectrum. For larger values of Ry, L changes 
slowly. Note that R C M is anti-correlated with L. As Ry 
increases, R C M grows sharply. But for Ry « 1000 (Pm ~ 
0.3) the threshold in R C M settles around 300, and the 
system seems to reach an asymptotic regime. 

The same quantities are shown in Fig. Q for the 
Roberts forcing (see also Ref. i)- In this case, the 
external forcing is an exact solution of the Euler’s equa¬ 
tion, and for Ry < 100 the flow is hydrodynamically 
stable. Between Ry « 100 and Ry « 1000 the flow 
displays hydrodynamic oscillations, and for Ry > 1000 
the flow reaches a hydrodynamic turbulent regime with a 
Kolmogorov’s power law in the inertial range. As a result 
of these instabilities, the characteristic scale of the flow 
changes twice sharply. In the laminar regime L ~ 27 t // cr . 
Then L decreases and remains approximately constant in 
the oscillatory regime (400 < Ry < 1000), and finally has 
a second sharp decrease {Ry « 1000) as turbulence de¬ 
velops. The threshold for dynamo action shows again an 
anti-correlation with L , and two sharp increases in R C M 
as a function of Ry are observed as L changes. In the 
turbulent regime {Ry « 1000), an asymptotic value of 
R c m seems again to be reached. In this case, only DNS 



R v 

FIG. 2: Integral lengthscale of the flow as a function of Rv 
(above), and R C M as a function of Rv (below) for Roberts 
forcing. 



Rv 

FIG. 3: Integral lengthscale of the flow as a function of Rv 
(above), and R C M as a function of Rv (below) for ABC forcing. 
Crosses connected with solid lines are obtained from DNS, 
while crosses connected with dotted lines are obtained using 
the LAMHD equations. 

was used. 

Figure [3] shows the threshold and integral scale as a 
function of Ry for ABC forcing. As mentioned in the 
previous section, this is a helical flow with energy injec¬ 
tion at intermediate scales, and as a result large scale dy¬ 
namo action can take place. As for Taylor-Green forcing, 
crosses connected with solid lines are obtained from DNS, 
while crosses connected with dotted lines are obtained 
from LAMHD simulations. There is an anti-correlation 
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FIG. 4: Volume render of enstrophy density (above) and 
square current density (below) in a small region of a simu¬ 
lation with Pm = 5 x 10 2 and ABC forcing. 


between L and Ry, but it is less clear than in the pre¬ 
vious flows. In the laminar case ( Ry ~ 10) the integral 
scale is close to 2tt/ kABC, while in the turbulent regime 
(Ry « 1000) R c m is larger and L is smaller. However, 
there is an intermediate range where L grows together 
with R c m - The increase of L in this intermediate regime 
is associated with the development of chaotic trajectories 
connecting several ABC cells as the flow bifurcates. The 
hydrodynamic bifurcations and subsequent generation of 
small scales (breaking the infinite time correlation of the 
laminar flow) also gives rise to an enhanced diffusivity 
that in turns increases R C M - 

When compared with the previous flows, there are two 
striking differences in the amplitude of R C M as a function 
of Ry. First, the thresholds at high Ry are reduced by 
one order of magnitude. Second, there is only a small de¬ 
pendence of R c m with Ry. As the flow destabilizes and 
develops turbulence {Ry « 100) the threshold changes 
from R c m ks 10 to R C M « 20, while in the previously dis¬ 
cussed flows the increase is by a factor of 10. After that, 
R c m drops and finally stabilizes. Note that the smallest 
magnetic Prandtl number of Pm « 5 x 10~ 3 was reached 
for this flow. 

Figure |1] shows renders of the enstrophy density and 
mean square current in a small box in a DNS using ABC 
forcing and Pm = 5 x 10 -2 . The scale separation be¬ 
tween the kinetic and magnetic diffusion scales is evident. 
While the velocity field displays thin and elongated vor¬ 
tex tubes, the structures in the magnetic field are thicker. 



FIG. 5: Spectrum of kinetic energy (thick line) and magnetic 
energy (thin line) in a 1024 3 simulation with Taylor-Green 
forcing (Pm =0.1 and Rm = 400). The spectrum of mag¬ 
netic energy has been multiplied by 10. Kolmogorov’s and 
Kazanstsev’s power laws are shown as a reference. 


In simulations without large scale flows El El no 
asymptotic regime for low Pm has been found so far. 
As Ry is increased and the system develops turbulent 
fluctuations, R C M grows to values much larger than the 
ones obtained here. From the numerical results, it is 
not clear whether an asymptotic regime will be reached 
or R c m will continue growing. However, theoretical argu¬ 
ments [3lj suggest that R C M should be constant for values 
of Ry large enough. From the comparison with the cases 
studied here and the anti-correlation found between L 
and R C M , it is clear that the presence of a large scale flow 
(helical or not) plays a crucial role in the development 
of the asymptotic regime at relatively small values of Ry 
and Rm- 


IV. SCALE INTERACTIONS 


The results presented in the previous section raise the 
question of what are the effects of turbulence, beyond 
increasing the value of R C M . Does dynamo action only 
take place at large scales when Pm < 1, or can turbulent 
fluctuations also amplify small scale magnetic fields? By 
small scales, we refer here to scales smaller than the flow 
integral scale but larger than the magnetic diffusion scale. 

Figure 0 shows the kinetic and magnetic energy spec¬ 
trum in the kinematic regime of a 1024 3 DNS using 
Taylor-Green forcing {Pm = 0.1 and Rm = 400). The 
magnetic energy spectrum peaks at scales smaller than 
the flow integral scale, and at intermediate scales the 
spectrum is approximately flat. At large scales, a slope 
compatible with Kazantsev’s spectrum is observed, 
although it should be noted that the hypothesis used by 
Kazantsev do not apply to this case. In Ref. Q it was 
shown that the peak in the magnetic energy spectrum 
moves to small scales as Ry is increased. 

The role played by the different scales in the flow is fur¬ 
ther clarified by examining the energy transfer functions. 
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FIG. 6: Transfer functions — Tl[P) (solid line) and Tjv/(fc) 
(dashed line) in the 1024 3 simulation with Taylor-Green forc¬ 
ing. 

From Eqs. CD and 0, the transfers 

mk)=J%- (r^B) k ^ k , (ii) 
T M (k) = J B£ ■ V x (^TB) k dfi k , (12) 

can be defined (see e.g. Ref. @). The asterisk de¬ 
notes complex conjugate, the hat Fourier transform, 
the subindex denotes the amplitude of the mode with 
wavevector k, and is the element of surface of the 
sphere of radius k in Fourier space. When negative, 
Ti(fc) represents energy given by the velocity field at the 
wavenumber k to the magnetic field at all scales. Positive 
Tm(^) represents energy received by the magnetic field 
at wavenumber k from the velocity field at all scales, and 
from the direct cascade of energy. The detailed shell-to- 
shell energy transfer in the kinematic dynamo has been 
studied in Ref. Jt>| . 

Figure |C| shows the transfer functions for the 1024 3 
DNS with Taylor-Green forcing. The magnetic field re¬ 
ceives energy from the velocity field in all scales in the 
inertial range (note that Tj, has constant amplitude from 
k ss 3 up to k « 40). As a result, for this forcing both 
the large scale flow and the turbulent fluctuations give 
energy to the magnetic field. On the other hand, Tm(/c) 
peaks at k ss 60. This suggests that the dynamo pro¬ 
cess is non-local in Fourier space. Other indications of 
non-local interactions from simulations were obtained in 
Ref. 1341 . and in simulations of dynamo action at large 
Pm I , although direct verification of the nonlocality in 
MHD turbulence was not obtained until recently :3a. 1371]. 

Similar results were obtained for Roberts forcing Q • In 
that case, a magnetic energy spectrum peaking at small 
scales and an approximately flat T b ( k ) transfer function 
at small scales were observed. 

The results for ABC forcing are different. Figure [3 
shows the kinetic and magnetic energy spectrum in a 512 3 
DNS with Pm = 1-5 x 10 -2 and Rm = 40. In this case, 



FIG. 7: Spectrum of kinetic energy (thick line) and mag¬ 
netic energy (thin line) in a 512 3 simulation with ABC forcing 
[Pm = 1.5 x 1CT 2 and Rm = 40). The spectrum of magnetic 
energy has been multiplied by 10. 


the magnetic energy spectrum peaks at k = 1 , and no 
peak is observed at scales smaller than the integral scale 
of the flow. Moreover, the magnetic energy spectrum 
drops fast. The transfer T^k) (Fig. 0 shows that most 
of the amplification is done by the large scale flow at 
k = 3, and the turbulent fluctuations seem to give only 
a small contribution to the dynamo. Tm{P) shows that 
most of the energy is received by the magnetic field at 
large scales, and it drops as T^[k) at small scales. 

The reason for this behavior can be understood as fol¬ 
lows. For ABC forcing, R C M is of the order of a few tens, 
even in the fully turbulent regime. As a result, a dynamo 
simulation just above threshold has small Rm , and small 
magnetic scales are damped so fast that no inertial range 
can develop in the magnetic energy spectrum. As a re¬ 
sult, the magnetic field is only amplified by the large scale 
flow and grows at larger scales. Since there is no dynamo 
source at scales smaller than the flow integral scale, the 
small scales are only fed by the direct cascade of energy 
and the magnetic energy spectrum cannot peak at these 
scales. 

However, if a simulation using ABC forcing is done at 
larger values of Rm, dynamo amplification by the small 
scales is recovered and similar features than for Taylor- 
Green forcing are obtained [2^,0]. This situation should 
be the one that prevails in astrophysics. 


V. CONCLUSIONS 

We reviewed results of dynamo action at low magnetic 
Prandtl number liQ for several mechanical forcing 
functions, including helical and non-helical flows, as well 
as small and large scale dynamo action. For all cases 
where a large scale flow is present, a similar behavior is 
obtained in the threshold for dynamo action: as Ry is 
increased (or Pm decreased), the value of R C M increases 
sharply as turbulence develops and finally reaches an 
asymptotic regime independent of the value of Pm- The 
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FIG. 8: Transfer functions —TL(fc) (solid line) and Tjv/(fc) 
(dashed line) in the 512 3 simulation with ABC forcing. 

large scale flow plays an important role in the establish¬ 
ment of the asymptotic behavior, as shown by the anti¬ 
correlation between the characteristic length of the flow 
and R c m - As turbulence develops, the laminar flow cre¬ 
ates small scales through hydrodynamic instabilities, and 
the large scale laminar flow loses its infinite correlation 
time. Then, both the large and small scale velocity fields 
amplify the magnetic field, giving rise to the asymptotic 
regime. 

New results from simulations with ABC forcing present 
a distinctive behavior. Having the flow maximum kinetic 
helicity and permitting large scale dynamo action, the 


critical magnetic Reynolds is smaller than for the other 
two flows by one order of magnitude. Also, only a twofold 
increase in R C M is observed as turbulence develops (in 
contrast to a tenfold increase for the other flows). As a 
result, dynamo simulations close to the threshold do not 
show small scale amplification (down to Pm = 5 x 10 3 ), 
and only the large scale flow is responsible for the large 
scale dynamo action. However, as the value of Rm is in¬ 
creased, as is expected for astrophysical and geophysical 
flows, small scale amplification in the kinematic regime 
is recovered. 

Although the conditions in our simulations are ideal¬ 
ized, we believe the existence of an asymptotic regime 
for Pm < 1 has profound implications for laboratory ex¬ 
periments and modeling of astrophysical and geophysical 
dynamos. In most of these cases, a large scale flow (such 
as differential rotation) and turbulent fluctuations are 
known to be present. 
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